Oh! 



New Ekpyrotic Cosmology 



o 
o 

(N 

> 
O 

^, 

^ , ^ Perimeter Institute for Theoretical Physics 

C\l ■ Waterloo, Ontario, N2L 2Y5, Canada 



Evgeny I. Buchbinder^, Justin Khoury^, Burt A. Ovrut^ 



^ ■ Department of Physics 

! The University of Pennsylvania 

Philadelphia, PA 19104-6395, USA 

(N ■ 
O 

o 



Abstract 



In this paper, we present a new scenario of the early Universe that contains a pre big 
bang Ekpyrotic phase. By combining this with a ghost condensate, the theory explicitly 
violates the null energy condition without developing any ghost-like instabilities. Thus 
^ , the contracting universe goes through a non-singular bounce and evolves smoothly 

^ . into the expanding post big bang phase. The curvature perturbation acquires a scale- 

invariant spectrum well before the bounce in this scenario. It is sourced by the scale- 
invariant entropy perturbation engendered by two ekpyrotic scalar fields, a mechanism 
recently proposed by Lehners et al. Since the background geometry is non-singular 
at all times, the curvature perturbation remains nearly constant on super horizon 
scales. It emerges from the bounce unscathed and imprints a scale-invariant spectrum 
of density fluctuations in the matter-radiation fluid at the onset of the hot big bang 
phase. The ekpyrotic potential can be chosen so that the spectrum has a "red" tilt, in 
accordance with the recent data from WMAP. As in the original Ekpyrotic scenario, 
the model predicts a negligible gravity wave signal on all observable scales. As such 
"New Ekpyrotic Cosmology" provides a consistent and distinguishable alternative to 
inflation to account for the origin of the seeds of large scale structure. 



1 Introduction 



Over the past decade observations of the microwave background temperature anisotropy 
have revealed that the large scale structure in our universe originates from primordial per- 
turbations that are nearly scale- invariant, adiabatic and gaussian. Since these coincide with 
the predictions of the simplest inflationary models, this is widely regarded as evidence for 
inflation. However this does not constitutes a proof, and it is prudent to keep in mind that 
the seeds for structure formation could originate from a different mechanism. Ultimately our 
faith in inflation must rely on the absence of a compelling alternative paradigm. 

In this paper, we present a fully consistent and complete scenario of early universe cos- 
mology which produces a nearly scale-invariant spectrum of density fluctuations without 
invoking a period of accelerated expansion. The model is strongly inspired by and borrows 
key ingredients from ekpyrotic [H [21 E] and cyclic [1] cosmology, in particular the idea that 
density perturbations are generated during a slow contracting phase prior to the big bang. 
A key difference, however, is that the cosmological evolution is now completely non-singular 
and avoids any big crunch singularity. Moreover, perturbations remain in the linear regime 
throughout. Thus, entirely at the level of a 4d effective theory, we are able to produce a 
scale-invariant spectrum of density fluctuations, track its evolution through the reversal from 
contraction to expansion, and its transmission to the matter-radiation fluid at the onset of 
the hot expanding phase. The idea of a contracting phase prior to the big bang originated 
in pre-big bang cosmology. Density perturbations and the evolution from the contracting to 
the expanding phase have been discussed in this context. See [5] for a review. 

The key ingredient in ekpryosis is a scalar field rolling down a negative, nearly exponential 
potential. If the coefficient in the exponent is sufficiently large, then, as the scalar field 
rolls down the potential, its fluctuations acquire a nearly scale-invariant spectrum spanning 
60 e-foldings in comoving wavenumber [H E]. Since the potential energy is negative, this 
occurs in a slowly contracting Friedmann-Lemaitre-Robertson- Walker (FLRW) geometry 
and, hence, with a rapidly decreasing Hubble radius. In some sense this ekpyrotic phase 
is "dual" to inflationary cosmology |6l [7j, which, in contrast, is a period of exponential 
expansion with nearly constant Hubble radius. The difference in dynamics, however, leads 
to a distinguishing observational prediction: the inflationary gravitational wave spectrum 
is nearly scale- invariant, whereas that of ekpyrosis is not [HE]- In ekpyrotic theories, the 
gravitational wave spectrum is strongly blue and, hence, the amplitude is exponentially 
suppressed on all observable scales [S]. 
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An important obstruction facing ekpyrotic theory is how to "bounce" from the contracting 
phase to an expanding phase, which requires a violation of the null energy condition (NEC). 
This is no small feat, however, since non-singular theories that violate the NEC generally 
suffer from violent instabilities, such as ghosts or tachyons of arbitrarily large mass |9]. In 
the Big Bang/Big Crunch ekpyrotic scenario of |2l [3], as well as in the cyclic model [1], 
one therefore allows the FLRW space-time to crunch, invoking stringy effects to effectively 
violate the NEC and generate a smooth bounce. This is motivated by the relative mildness 
of the singularity [2j. (Understood, for example, as a collision of end-of-the- world branes in 
heterotic M-theory [lOl [11] , only the fifth dimension shrinks to zero size while the 3 large 
dimensions remain finite.) Despite considerable effort in string theory [12] and compelling 
physical arguments [13J, there is still no proof that a bounce is possible. 

Recently, however, NEC-violating solutions [H] have been derived in the context of ghost 
condensation [151 IIS] • Since these models involve higher-derivative kinetic terms, they evade 
the assumptions of [9] and therefore yield ghost-free solutions. 

In this paper, we show explicitly how the NEC-violating ghost condensate can be merged 
consistently with the preceding ekpyrotic phase to generate a non-singular bouncing cos- 
mology. This merger involves many subtleties. For instance, since the energy density of an 
NEC-violating fluid grows in an expanding universe, whereas everything else redshifts, it 
quickly comes to dominate the universe. But in a contracting universe precisely the oppo- 
site happens: the NEC-violating component goes to zero while any amount of radiation or 
scalar kinetic energy blueshifts. Thus achieving a bounce hinges on an efficient transfer of 
energy from the ekpyrotic scalar, which dominates during the ekpyrotic phase, to the ghost 
condensate, which achieves the bounce. To realize this we therefore propose that they are, 
in fact, one and the same field. The higher-derivative kinetic function is chosen to be nearly 
canonical during the ekpyrotic phase, while higher derivatives become relevant only as we 
approach the bounce. This translates into consistency relations between the kinetic function 
and the scalar potential, which we derive explicitly. We argue that these conditions are 
realized for a wide class of kinetic functions and scalar potentials. 

Armed with a non-singular cosmological evolution, we can address the propagation of 
density perturbations through the bounce, an issue that has stirred a lot of controversy over 
the past few years [17J. Although the fluctuations in the scalar field and the gauge-invariant 
Newtonian potential are scale-invariant, it turns out that their growing mode precisely can- 
cels out of C — the curvature perturbation on uniform- density hypersurfaces [I8] . The latter 
is a useful variable to track since it is conserved on super-horizon scales, in the absence 
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of entropy perturbations. Thus, as long as 4d general relativity is valid, ( remains con- 
stant independent of the physics of the bounce. In the context of the Big Bang/Big Crunch 
ekpyrotic or cyclic models, however, 4d gravity must break down near the singularity and 
therefore the issue of what happens to ( remains unresolved. Some maintain that the likely 
matching condition is for ( to be continuous, most convincingly [19], while others argue that 
higher-dimensional effects near the singularity could lead to mode-mixing and endow ( with 
a scale-invariant contribution [20l [211 122] . Nevertheless, one must also deal with the fact 
that this mode diverges logarithmically near the singularity, resulting in a breakdown of 
perturbation theory, although matching prescriptions based on analytic continuation have 
been proposed [201 123] • 

In our case, the evolution is completely non-singular and, hence, the outcome of the 
perturbations is unambiguous: ( is conserved. To generate a scale-invariant spectrum we 
rely instead on a recently proposed mechanism using entropy perturbations [1211 [25]. For 
earlier and closely related work, see f2E[ f27\ [25 1 If we have two scalar fields, each rolling 
down a steep, negative, and nearly exponential potential, then each of them acquires a scale- 
invariant spectrum of fluctuations, as described earlier. Moreover, the entropy perturbation 
— corresponding to the difference in the scalar fluctuations — is also scale-invariant [2^ [27] . 
By converting this entropy perturbation into the adiabatic mode {i.e., (), the curvature per- 
turbation thus becomes scale-invariant long before the bounce. In the inflationary context, 
the idea of using a spectator field to generate a scale-invariant spectrum of entropy perturba- 
tions, later to be imprinted on the adiabatic mode, was also considered in the curvaton [30] 
and modulon scenarios [21]. 

In this paper, we propose a concrete realization of the entropy generation mechanism. 
We consider two ekpyrotic scalar fields, each with its own negative potential and its own 
higher-derivative kinetic function. As the fields roll down their respective potentials, the 
entropy perturbation acquires a scale-invariant spectrum. To convert this spectrum into the 
adiabatic mode, we assume that the ekpyrotic potentials are such that one field exits the 
ekpyrotic phase and enters the ghost condensation regime before the other. More precisely, 
the potential in the ekpyrotic phase is steep and negative, as mentioned earlier, while in the 
ghost condensate phase it must be positive. Thus this transition is marked by a sharp rise in 
the potential. If the potentials are such that one field hits the transition before the other, this 
will result in a sharp turn in the trajectory in field space, and in the process ( will acquire a 
scale-invariant piece from the entropy perturbation. We calculate the resulting ( explicitly 
and find an expression that closely resembles its inflationary counterpart. This confirms 
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earlier claims that the required level of tuning on the ekpyrotic potential is comparable to 
that in inflation. 

For exact exponential potentials, we flnd that the resulting spectral tilt is slightly "blue" , 
which at flrst sight is disturbing in light of the recent WMAP evidence for a small "red" 
tilt. However, we show that deviations from the pure exponential form can lead to a small 
red tilt. See also [21] for an independent derivation. This is shown explicitly in the case 
that the potentials for the two scalar flelds are identical, albeit not exactly exponentials. 
This simpliflcation is made only to facilitate the analysis, and we do not believe that the 
resulting red tilt hinges on it. In this limit we derive an exact evolution equation for the 
entropy perturbation, cast entirely in terms of the background equation of state, which closely 
resembles analogous results in inflation [32] and old (single-fleld) ekpyrotic theory [3^ . The 
resulting spectral tilt acquires a dependence on the degree of departure from pure exponential 
form, and, in particular, can be red. 

Coming back to the background evolution, one of the potential dangers with contracting 
universes is that the FLRW background is unstable to the onset of chaotic mixmaster be- 
havior [3H |35] . In [36] it was shown that mixmaster behavior is suppressed if the dominant 
energy component has equation of state if ^ 1, as is the case in the ekpyrotic contracting 
phase. In singular ekpyrotic theories, as the universe approaches the singularity the equation 
of state must eventually revert to w = 1, and mixmaster behavior can potentially resurface. 
A necessary condition to avoid a chaotic bounce is that the anisotropy and curvature com- 
ponents be exponentially suppressed at the onset of the w = 1 phase so that they remain 
subdominant all the way to the string or Planck scale, where 4d gravity breaks down anyway. 

Mixmaster issues also apply to our new scenario but are trivially satisfled since the 
bounce is non-singular. The w ^ 1 ekpyrotic phase exponentially suppresses anisotropy and 
curvature components. The latter grow again during the NEC-violating phase, but only by 
an insignificant amount since the bounce occurs within one e-fold of contraction. 

To summarize, our scenario is a complete and unambiguous template for early universe 
cosmology. It provides an alternative explanation for the origin of the seeds of structure for- 
mation in which perturbations are generated long before the bounce, all within a singularity- 
free and finite FLRW background. The bounce is entirely described at the level of a consistent 
effective theory, which is free of ghost-like instabilities or other pathologies. Mixmaster be- 
havior is trivially avoided. Perturbations remain in the linear regime, and their evolution can 
be tracked through the bounce. Since 4d gravity remains valid throughout, the curvature 
perturbation ( goes through unscathed. And since ( acquires a scale-invariant spectrum 
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long before the bounce, thanks to the entropy perturbation conversion mechanism, it is un- 
equivocally scale-invariant after the transition. The universe therefore emerges in a hot big 
bang phase, endowed with a super-horizon spectrum of scale-invariant fluctuations. While 
the approach presented here is closely related to and borrows key elements of the original 
ekpyrotic scenario [1], our new scenario resolves all of its short-comings, such as generat- 
ing a bounce. It also resolves the important issues of singularity-avoidance and the fate 
of the perturbations that occurs in [21 [3l Hj. For these reasons, we call it "New Ekpyrotic 
Cosmology" . 

As it stands the new ekpyrotic scenario is precisely that, a scenario. In principle, it might 
be implemented in various ways in different fundamental theories of particle physics. That 
being said, this scenario maintains the motivations of the original ekpyrotic model [I] , namely, 
as the cosmology associated with the singularity free collision of a bulk five-brane EB] 
in heterotic M-theory with the observable boundary wall. Indeed, we have computed the 
potential in such a theory and find that it can satisfy many of the constraints required 
in the new ekpyrotic model. This potential will be presented elsewhere ^39j. A further 
motivation is that a realistic matter spectrum can appear on the observable wall of such 
theories, see, e.g., [lO]. Hence, our scenario can potentially occur in a realistic context. 
More generally, we emphasize that, with the exception of the so-called no-scale theories, a 
supersymmetric minimum in the scalar space of any 4d = 1 supergravity theory, including 
effective low energy superstring theories, has a negative cosmological constant. This can 
remain true in the presence of flux and non-perturbative effects, even for a vacuum that 
breaks supersymmetry [Hj. Therefore, potentials consistent with ekpyrotic cosmology, i.e., 
where scalars roll down a steep negative potential energy, appear naturally in this context. 

A more pressing question is whether a ghost condensate can occur in string theory. It 
has been argued that this is impossible given the known analytic properties of the string S- 
matrix [l2]. See [13], however, for a recent attempt to overcome these difficulties. Without 
question higher derivative interactions do occur in string theories, but whether they are of 
the requisite form is an open issue. Ultimately a non-singular bounce might arise from very 
different physics, unrelated to string theory. In fact, the NEC could conceivably be smoothly 
violated by a different mechanism than ghost condensation. Be that as it may, many of our 
results would continue to apply to any such mechanism. The rationale for focusing here on 
the ghost condensate is to provide an explicit and complete realization of our new scenario. 

In Sec. [2] we review the essential concepts of the ekpyrotic phase, as a theory of a scalar 
field rolling down a negative potential. We show how the fluctuations in this field acquire 
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a scale-invariant spectrum, which unfortunately projects out of the curvature perturbation. 
In Sec. [3] we extend the discussion to two scalar fields and review the mechanism proposed 
by [Ml |25] to generate a scale-invariant spectrum for the entropy perturbation. In Sec. H] 
we derive an explicit expression for the spectral tilt and argue that it can be slightly red, 
in agreement with recent microwave background data. We show in Sec. [5] how this gets 
imprinted into the curvature perturbation by requiring that the two fields exit the ekpyrotic 
phase at different times. Next we turn to the bounce, first, in Sec. El discussing general 
requirements on its physics followed by a brief review of ghost condensation in Sec. [3 In Sec. [8] 
we show in detail that the ekpyrotic phase and ghost condensate can be merged successfully 
to generate a non-singular cosmological scenario, and derive consistency requirements on the 
kinetic function and scalar potential. Section [9] presents a short discussion of the physics of 
reheating in this scenario, while Sec. [TO] provides some concluding remarks. 

2 Review of Single-Field Ekpyrosis 

At the level of a 4d effective description, the basic ingredients of the simplest ekpyrotic 
scenario are essentially the same as in inflation, namely a scalar field rolling down some self- 
interaction potential V{(j)). A key difference, however, is that while inflation requires a fiat 
and positive potential, its ekpyrotic counterpart is steep and negative. This has a dramatic 
impact on the cosmological evolution. Instead of accelerated expansion, an ekpyrotic theory 
has slow contraction. Instead of an exponentially growing scale factor and nearly constant 
Hubble radius, corresponding to approximate de Sitter geometry, we now have a nearly 
constant scale factor and rapidly shrinking Hubble radius, corresponding to approximately 
fiat space. 

2.1 Ekpyrotic Potential 

A generic ekpyrotic potential, shown in Fig. [T|, consists qualitatively of three parts. In the 
region denoted by (a), the potential must be steep and negative. As the field rolls down this 
part of the potential, there is a scaling solution, which is an attractor, corresponding to very 
slow contraction. 

It is also during this phase that large-scale density fluctuations are generated as modes 
exit the horizon. As we will review later, in order for the spectrum to be nearly scale- 
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Figure 1: Generic shape of the ekpyrotic scalar potentiaL 



invariant, the potential must satisfy the "fast-roll" conditions 



e < 1 ; |?7| -Cl 



where 



e = M 



PI 



V 



(2.1) 



(2.2) 



are fast-roll parameters, in analogy with the standard slow-roll parameters in infiatioiij. 
Here = dV/dcp, and Mpi is the "reduced" Planck mass: Mpi = 2.4 x 10^*^ GeV. These 
conditions respectively require the potential to be steep and nearly exponential, and thus 
region (a) of the potential can be approximated by 



V{(j)) ^ — Voexp ( — 



(2.3) 



with p ^ 1 . 

Much as inflation ends when the flatness condition breaks down, here as well the scaling 
behavior terminates once equations ( 12. 2p are no longer satisfied. In order to avoid being 
left with a large negative vacuum energy at the end of the ekpyrotic phase, let us assume 
that the potential has a minimum and rises back up towards positive values, as shown in 
region (b) of Fig. 1. 



^Recall that the slow-roh parameters in inflation are einf = ('^) " "^^Pi i~v^\ 
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There is considerable freedom in specifying the shape of the potential further to the left of 
the minimum, that is, in region (c). Since it is in this region that the universe reverses from 
contraction to expansion, the detailed shape of the potential is dependent on the explicit 
mechanism producing the bounce. In this paper, we will use ghost condensation to generate 
a non-singular bounce where all instabilities are under control. As we will see in Sec. U\ this 
requires a flat and positive potential in region (c), as sketched in Fig.l. 

2.2 Scaling Solution 

The ekpyrotic phase occurs as the field rolls down the steep, negative exponential part of 
the potential; region (a) in Fig. [H In this paper, we take the background geometry to be 
a homogeneous, isotropic Friedmann-Lemaitre- Robertson- Walker (FLRW) space-time which 
is spatially flat. That is, 

ds'^ = -df + aitfdx^. (2.4) 

Furthermore, in this section, is chosen to have canonical kinetic energy and no higher 
derivative interactions. The equations of motion are then given as usual by the Friedmann 
equation, 

m^Ml, = + V(0) , (2.5) 
as well as the evolution equation for the scalar field: 

4> + 3H<i)= -V^. (2.6) 

For the potential (12.31) . it is easily seen that these equations allow for an exact scaling 
solution [31 [331 Hg 

a{t) ~ i-tf ■ H = ^; 



0(t)^V^Mp.log(^-^ ^,J°_3^^ tj, (2.7) 

where t is negative and increases towards zero. Hence the solution describes a slowly con- 
tracting universe. Since p <^ 1, we see that the kinetic and potential energy of the scalar are 
both large in magnitude, but nearly cancel each other to yield a small total energy density. 
In other words, <j)'^/2 > H^M^^, and \V\ > H^M^y. 

This scaling solution, moreover, has the desirable property that it is an attractor. Indeed, 
the scalar field has the equation of state of a very stiff fluid: 

V 2 

w = - = 1>1. (2.8) 

P 3p 



This means that its energy density behaves as ~ a~^/^. But since p <C 1, this energy 
component blueshifts much more rapidly than any other relevant contribution to the Fried- 
mann equation. Any curvature (~ matter (~ radiation (~ a~^), or anisotropy 

(~ a~^) component, for instance, quickly becomes subdominant to the scalar field energy 
density. Incidentally, this is precisely analogous to inflation where the scalar potential gives 
a nearly constant source in the Friedmann equation. Such a term therefore redshifts much 
more slowly than any of the aforementioned relevant contributions, which is why accelerated 
expansion is an attr actor. 

2.3 Density Perturbations — Newtonian Potential Analysis 

By far the most important feature of a negative exponential potential, however, is that 
it generates a scale-invariant spectrum of fluctuations in the scalar field, even when the 
gravitational interactions are turned off. Indeed, in the absence of gravity, (12.61) can be 
integrated trivially: 

= -V-2V , (2.9) 

where we have set the total energy to zero. For the pure exponential potential (12.31) . the 
solution is 

Now fluctuations in the scalar field satisfy the perturbed equation 

+ = 0, (2.11) 

which, using (I2.10p . can be written as 

50.+ 1)^0^ = 0. (2.12) 

This describes a harmonic oscillator with a time-dependent mass. Assuming standard Bunch- 
Davies boundary conditions [15], Scj)^ — > e^^^ j^J^ as k(—t) oo, the general solution for 
the mode functions is given by a Hankel function of the first kind: 6(j)k ~ \/^H'^j'^{—kt) . In 
the long- wavelength limit, k{—t) 0, this gives 



k-% (2.13) 



corresponding to a scale-invariant power spectrum. Scale invariance can be traced back to 
the factor of 2 coefficient in the time-dependent term in (I2.12p . Remarkably, from (12.101) we 
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see that this holds for any p. When we include gravitational interactions, however, we will 
find that scale invariance requires p <^ 1. 

Turning on gravity, the gauge-invariant variable that faithfully reproduces the above 
scale-invariant spectrum is the Newtonian potential $, the scalar metric fiuctuation in New- 
tonian gauge. It is convenient to do the analysis in terms of a related gauge-invariant variable, 
u = a^/(j)', where primes denote differentiation with respect to conformal time r. Its Fourier 
modes satisfy [3] 

< + (k' - TT-^TJ-j) u, = (2.14) 



whose solution is once again given by a Hankel function 

where z/ = (1 +j9)/2(l — p), and we have assumed Bunch-Davies initial conditions: Uk 
g-«fcT^^2/^)3/2^ The power spectrum on large scales is then 

eul ~ k-^r>/{i-P) ^ (2.16) 

which is indeed scale- invariant if p -C 1. The spectral index in this limit is ra^ — 1 ~ — 2p, 
corresponding to a slight red tilt. For more general potentials, the expression for in terms 
of the fast roll parameters (12. 2p closely resembles its slow-roll counterpart in inflation: 

ns-l = -A{e + r]) . (2.17) 

In particular, potentials with non-trivial ?7, characterizing deviations from pure exponential, 
can lead to a blue tilt. 

The nearly static nature of the geometry during this ekpyrotic phase is underscored by 
the fact that only a few e-folds of contraction are needed to generate 60 e-folds worth of 
perturbations. The comoving wavenumber for any given mode is related as usual to the 
moment of horizon crossing by ~ a\H\ ~ (— 1)^~^. Letting i and / denote respectively the 
initial and flnal time during which perturbations are generated, we therefore have 



ki \tf I \a. 



Since the spectral tilt is constrained observationally to be within 5% of exact scale- invariance, 
we have p ~ 1/40, and therefore 

In (^^^^3/2. (2.19) 

Thus, in contrast with inflation where the universe grows exponentially big as the relevant 
range of modes is generated, here the universe only shrinks by a factor of order unity. 
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2.4 Curvature Perturbation 



While u is all we need to specify the scalar metric perturbations, it is useful to keep track 
of another gauge-invariant variable C, — the curvature perturbation on uniform-density hy- 
persurfaces. This variable has the virtue of being constant on super-horizon scales, since 
there is no entropy perturbation in this single-field case. Thus, the large-scale spectrum of C 
calculated during the scaling phase automically agrees with the spectrum at horizon re-entry 
after reheating. 

Unfortunately, it turns out that the scale-invariant growing mode of $ calculated above 
precisely projects out of C- Since the two are related by 



it is easy to check using fl2.15p and the background solution that the Newtonian potential 
precisely behaves as $ ~ a' /a? on large scales, leaving C, with an unnacceptably strong blue 
tilt. This cancelation is at the core of all the controversy about the ekpyrotic spectrum of 
perturbations. 

While the derivation is rigorous, the conclusion is not without caveats. A crucial as- 
sumption is that the physics of the bounce from contraction to expansion all lies within the 
regime of validity of 4d effective theory, since the statement about C, remaining constant 
relies on the equations of 4d gravity and matter. The possibility that higher-dimensional 
effects relevant at the bounce could imprint a scale-invariant contribution to C, has been the 
subject of ample literature [20l [22] . 

In this work, however, we shall exploit an alternative possibility relying on two scalar 
fields and based on an observation by N. Turok [25] and work of Lehners, McFadden, Stein- 
hardt and Turok |2l] . This mechanism was suggested by [26] and exploited in [271 [2H1 [29] . 
Namely, with two scalar fields each rolling down their steep exponential potential, the entropy 
perturbation acquires a scale-invariant spectrum on large scales. If this can subsequently be 
imprinted onto then C, is endowed with a scale-invariant spectrum well before the bounce. 
After reviewing the mechanism of [231 [25] , we will show how the entropy perturbation im- 
prints a scale-invariant contribution to C, if one of the two fields reaches the minimum of its 
potential (region b) of Fig. [T]) before the other. Since this field loses much kinetic energy in 
the process, this corresponds to a sharp turn in the trajectory in field space, and triggers a 
scale- invariant jump in C,. Then, by the same argument as above, whatever bounce physics 
comes into play to reverse contraction to expansion, as long as everything remains 4d and 
perturbative, the scale-invariant C, will go through unscathed. 




(2.20) 
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3 Two-Field Ekpyrosis 



Remarkably, the scaling solution described above generalizes to two fields [27] . Consider 
two scalars and ip, each with a potential of the form shown in Fig. [H In region (a) we can 
approximate the scalar potential as a sum of exponentials, generally with different powers 
p and q <^ 1: 

Vi^M = -Voexp (-^1^) - f/„exp {-^^^) ■ (3.1) 

Furthermore, we continue to assume that each scalar field has canonical kinetic energy and 
that there are no higher derivative interactions. The dynamics of and ip are therefore 
governed by, respectively, 

4> + 3H^= -V^ ; 

'ip + 3H^ = -V^ . (3.2) 

When we come to perturbations, it will be useful to describe the field trajectory in (0, ip) 
space in terms of the adiabatic field a, given by the geometrical relation 

CT = cos 6^0 + sin 6''0 , (3.3) 

where 

tan9 = i. (3.4) 



Combining (13.21) . it is easily seen that a satisfies 

a + 3Ha = -V„ , (3.5) 
where the slope of the potential along the field trajectory is just 

= cos^l^^ + sin^l/^. (3.6) 

In the exponential regime, there is an exact scaling solution to the Friedmann and scalar 
equations which generalizes (12. 7p to 

p + q 



a{t) ~ (-t)*'+« ; H 



t 



^<<) = V2.Mp,.og|-J^^^-^^A_,|. (3.7) 
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Unlike the single-field case, however, the above solution is not an attractor because of an 
instability in the direction orthogonal to the field trajectory. It is precisely this instability 
which is exploited in the next subsection to amplify entropy perturbations. Numerical anal- 
ysis reveals that the instability grows fastest along the direction of the steepest exponential. 
Suppose p < q so that the potential for (p is steeper than that for ip. Then, the instability 
brings ip to a. halt, while the solution converges to the single-field scaling solution fl2.7p along 
the direction. Nevertheless, this is not of great concern since one only needs to be on the 
scaling solution for a few e-folds of contraction to generate the observable range of modes, 
as (I2.19P demonstrates. 

3.1 Generation of Scale- Invariant 6s during the Scaling Regime 

In this subsection, we discuss the generation of a scale-invariant spectrum of entropy per- 
turbations for the above scaling solution [2H |27] . See [16] for a nice review of entropy and 
adiabatic perturbations in multi-field models. A convenient gauge for the study of multi- 
field perturbations is the so-called spatially-flat or Mukhanov- Sasaki gauge [17], in which 
the perturbations 6(f) and 64' in the two scalar fields satisfy the coupled equations 

^ ,HS, . . {k« - ^1 (p) } - {p'i') 

St + 3HH. + + {k** - {^^) ] S„ - (^h) S,^0. 

Thus, even though the scalar fields are uncoupled at the level of the potential, in the sense 
that V^if) = 0, their perturbations become intertwined through gravity. 

These fluctuations can be decomposed as a component along the field trajectory, 6a, 
which is the adiabatic perturbation, and a perturbation orthogonal to the trajectory, 6s, 
which is the entropy or isocurvature perturbation: 

6a = cos 9 6(f) + sin 9 6ip 

6s = - sm9 6(1) + cos 9 6ip. (3.8) 

Combining the {6tl),6(f>) equations, and using the standard energy and momentum constraints, 
the entropy perturbation can be shown to satisfy [16] 

6s + 3H6s +('^ + V,s + S9A 6s = 4M|,-^^ , (3.9) 

where \1/ is the curvature perturbation in Newtonian gauge, while Vss = cos^ 9 V^^ + 
sin^ 9 V^(f) is the curvature of the potential orthogonal to the field trajectory. 
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During the ekpyrotic scaling phase, described by (13. 7p . both and ip scale as 1/t and, 
thus, 9 is constant: 

tane=J-. (3.10) 
V P 

Furthermore, for the scaling potential (13.11) we have 

Vss = ~{l-3{p + q)) . (3.11) 
Plugging all of this into (13. 9p , we find 

6s + 3H6s + l — --{l-3{p + q))\6s = 0. (3.12) 

This can be simplified further by introducing a rescaled variable v = aSs and rewriting in 
terms of conformal time, dr = dt/a dtf^^^'^^: 

v" + ev - ^ (^l-^{p + q)^v ^0, (3.13) 

where we have dropped higher order terms in p, q. Here r is assumed negative and increas- 
ing towards zero. Since p and q are small, the coefficient of the time-dependent term is 
approximately equal to 2, indicating a nearly scale-invariant spectrum. Indeed, the solution 
to (I3.13P with standard Bunch-Davies initial conditions is, up to a phase, 

i=y|v^if«(-A:r), (3.14) 
with 



3 / 4, , 3 , , , , 

^=2V^"3^^ + ^^'^2"^^ + ^^- ^^-^^^ 
On large scales, k -C aH, the amplitude tends to w ~ k^", from which we can read off 
the spectral index: 

-1 ^2(j9 + g) =4e, (3.16) 

where we have generalized the fast-roll parameter e given in (12.21) to the two-field case in the 
obvious way: 

e = Mpi2(-^) <1. (3.17) 



This parameter measures the steepness of the potential along the field trajectory. Thus, for 
pure exponential potentials, the spectrum is slightly blue. More general potentials allow for 
red spectra as well, as we will show explicitly in the next section. 
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4 Spectral Tilt for General Potentials 



The small blue tilt for pure exponential potentials is a slightly disconcerting result in light of 
the recent 3- year WMAP data which favors a red tilt. However, we expect — and will show — 
that this unambiguous blueness is an artifact of the pure exponential form. This intuition 
is supported by the single-field result: a non-zero value for rj in fl2.17p can change the sign 
of the tilt. We will see that the two-field story is nearly identical. 

More general potentials immediately imply departure from the scaling behavior studied 
earlier and thus generically require numerical analysis. For simplicity we focus here on the 
case where both scalar fields have identical potential V in region (a) of Fig. [1) 

V^(0,V^) = V(0) + V(^). (4.1) 

This approach was originally mentioned in |i25|, and some of the results below overlap 
with Of course the potentials should only be identical in region (a), not globally, 

for otherwise the entropy perturbation will never get converted to the adiabatic mode. We 
can assume, for instance, that the minimum in region (b) occurs at different field values. 

If the fields start out with identical initial conditions, then their time-evolution will also 
be the same, and thus 6 = tt/A. It follows that the adiabatic field introduced in (13. 3p satisfies 
a = (0 + ■?/')/ V2. Similarly for the slope of the potential along the field trajectory, given 
by (13. 6p : Va = (V<^ + V^) /V^. Thus the Friedmann and "if" equations can be expressed 
as 

3H^Mi, = h'' + V; 

HMl, = -]^a\ (4.2) 

Combined with the a equation of motion (13.51) . this rewriting makes manifest the virtue of 
having identical potentials for the two scalars — all background equations of motion reduce 
effectively to that of a single effective scalar field a. 

Moving on to perturbations, the evolution equation ( 13. 9p for the entropy mode greatly 
simplifies because ^ = 0. Moreover, we can rewrite the time-dependent potential term, V^ss-, 
in terms of cr-derivatives: V^ss = cos^6'V^^ + sin^6'V0(/, = Vaa- Our starting point for the 
generation of Ss is therefore 

6s + 3H5s + + Ss = 0. (4.3) 
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Since the background no longer follows some exact scaling solution, the way to proceed 
is to recast every term in the equation solely in terms of the background equation of state 
parameter, 

e-2(l + -) = -^ = -^, (4.4) 

and its derivatives. Here = Ina is the number of e-folds, as usual. (In the limit of pure 
exponential potentials, e is of course constant and related to the e parameter in fl3.17p by 

6=1/26.) 

To illustrate the method, let us combine equations (14.21) to obtain a relation for the 
potential: 

V = H'^Ml.^ (3 - e) . (4.5) 

Then, from the second of (14. 2 p and the definition of e, we immediately have da/dN = 
\/2iMpi, which allows rewriting the derivative of the potential as 

V. = -H^M,,V2-e (^3 - e + . (4.6) 

Similarly, we have 

Vaa r)-2 p.- 5 (ie 3 dine 1 f dhie\^ Id'^lne 

ip^ + 6e + -— - - - \^-^ j - --^ . (4.7) 

Inspired by analogous calculations in the inflationary context [32J, it is convenient to 
rewrite (14. 3p in terms of a dimensionless time variable 

e — 1 aH 



and, as we did in Sec. 13. ![ rescale the perturbation variable to 

v = a5s. (4.9) 
After some algebra, the perturbation equation (14.31) takes the exact final form 

'rfln(e-l)V d^\n{e-l) 



e 



1 rfln(e-l)y ^d^v 



1 dN J dx'^ (e-iy 



dN / dN^ 



dv 2 

X- — \- X V 

dx 



if . 5 de 3dlne ifdlneV Id^lnel 

+ (F3Tp|-2-2^" + 7.-+-^--— --(^) = (4.10) 

As advocated, this depends only the background equation of state parameter and its 
time derivatives. To simplify things, we use the fact e is large during the ekpyrotic phase, 
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corresponding to w S> 1. Moreover, we assume it is a slowly- varying function of over the 
observable range of modes. To leading order, it follows from fl4.5p and fl4.6l) that 

(4.11) 

where e was defined in (13 .170 . Similarly, it is easy to show that 

dine , , , . 

^ = «,. (4.12) 

where we have generalized the rj parameter of (12.21) to the two-field case: 

^-1-^- (4.13) 

In the approximation that e and r] are small and nearly constant, (14.101) collapses to the 
simple form 

^'^4 + - 2 (1 - 3(e - r^)) = , (4.14) 

ctx^ 1 — or] 

where the coefficient of x'^ can be reabsorbed by a constant rescaling of the time variable. 
This final form can be recast as a Bessel equation, as before, and the corresponding spectral 
tilt can be read off immediately 

- 1 ^ 4 (e - r]) , (4.15) 

As a check, this agrees with (13.161) in the limit of pure exponential potentials. 

More importantly, we see that a red tilt is possible if t] is positive and dominates over 
the e term. This mild condition is satisfied by a wide class of potentials. For instance, any 
potential of the form V^(0) ~ exp(— 0"), with n > 2, gives a red tilt at large 0. 



5 Converting Entropy to the Adiabatic Mode 

Next we turn to the conversion of the scale-invariant entropy mode to the curvature pertur- 
bation. Indeed, the upshot is that the entropy perturbation sources C on large scales: 

C^-2H-5s. (5.1) 
a 

If the field trajectory is a straight line, so that 6 remains constant, then ( remains constant 
as well. If, however, the field trajectory has curvature or makes a sharp turn, then ( will 
change on large scales. In the context of our scenario, we present an explicit mechanism for 
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converting the nearly scale-invariant entropy perturbation into the adiabatic mode, which 
naturally exploits the desired shape of the potential shown in Fig. [TJ 

During the scaling solution, region (a) of the potential, 6s acquires a scale-invariant 
spectrum on large scales, as described earlier. Since this scaling solution corresponds to 
6 = const., however, it follows from (15. ip that ( is not scale- invariant, as in the single-field 
ekpyrotic model. Eventually the potentials for (p and ip must eventually depart from pure 
exponential by hitting a minimum and growing to positive values, corresponding to region 
b). The scaling solution abruptly ends. In general, one of the two fields, say ip, will hit 
the minimum before the other. As ip climbs up the steep hill, its kinetic energy decreases 
tremendously, while that of (p stays nearly the same. Hence, in field space this corresponds 
to a sharp turn in the trajectory, from 6 ^ arctan a/ q/p to 6 ^ 0. This jump in 6 implies a 
jump in ( jumps as well, which thereby acquires a scale-invariant spectrum from the entropy 
perturbation. 

If the rise in the ip potential is sufficiently steep, then the change in 6 occurs almost in- 
stantaneously compared to a Hubble time. This rapid-transition approximation is consistent 
with the fact that, in the scaling regime, we have ip ^ H since g <^ 1, as seen from (13. 7p . 
Thus, if the rise in the potential is indeed sufficiently steep, the field climbs on the plateau 
of region c) in a very short time compared to the Hubble time. Thus we can rewrite (15. ip as 



where we have denoted the time of the transition by tj. We can argue that each factor in this 
expression, save for the delta-function term, is approximately constant during the transition. 
This is obviously the case for H, by assumption. To see that the same is true of a, let us 
look back at (13. 5p . As the field trajectory makes a sharp turn, the slope of the potential 
Vcr will generically be discontinuous, i.e., have a theta-function jump. If & were to jump as 
well, however, the a term in (13.51) would generate a delta-function in the equation of motion, 
which is inconsistent. Hence a must be continuous, and we can therefore substitute its value 
of \/2{p + q)/ U just before the transition. Finally, 6s does change, but by at most a factor of 
order unity. The proof requires some steps, which we will provide at the end of the section. 
For the moment let us treat it as essentially constant. 

Making use of this rapid-transition approximation, we can integrate (15.21) and find that 
the curvature perturbation ( inherits a nearly scale-invariant spectrum from 6s: 






(5.3) 
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with spectral tilt given by (13.161) . Ignoring the change of order unity of during the transition, 

we can substitute the large-scale 6s by taking the hmit — > of (13.141) . neglecting for 
simplicity the small departure from scale-invariance: 

k^/^Ss ^ -J^ ^ ~^ . (5.4) 

v^a(-r) V2ip + q) 

Substituting this, as well as the scaling solution a = \/2HM-p\/ + g, we obtain 

~ arctan^ { , (5.5) 



~ arctan \ ^ — 

where e is the fast-roll parameter defined in (13.171) . Remarkably, up to the trigonometric 
factor of order unity, this is nearly identical to the corresponding expression in slow-roll 
inflation: 



where einf = M|i(V0/V)^/2 is the standard slow-roll parameter. This confirms earlier expec- 
tations that the level of tuning on the ekpyrotic potential from the amplitude of density 
perturbations is comparable to its inflationary counterpart. 

It remains to prove that the entropy perturbation 6s does not change dramatically dur- 
ing this process. Since the modes of interest are already far outside the horizon, we can 
study (13. 9p in the limit A; — > 0: 

6s + "iHSs + {v^ss + 3^^) 5s = . (5.7) 

For concreteness we consider the case where hits the minimum of its potential, climbs up a 
steep hill and hits a fiat plateau, as described above. Since loses most of its kinetic energy 
in the process, we have 9^0. Moreover, since the potential becomes fiat. Hence, 

V^ss — > as well. Finally, since the transition is assumed nearly instantaneous on a Hubble 
time, we can safely ignore the Hubble damping term. Thus (15. 7p reduces to: 

6s ^ -3e^6s . (5.8) 

Here the approximation 9 6{t) made earlier is too drastic since one would have to make 
sense of 6'^{t). Thus we have to be a more precise about 9{t) during the transition in order 
to solve for 6s. To proceed analytically, let us model the transition by a linear extrapolation 
between the initial and final angles, 9i and respectively: 

9{t)=9,+ ^^^{9f-9,). (5.9) 
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Thus the transition begins at t = tj and ends at t = tj + T, with HT -C 1. In the end we 
will choose 9i ~ arctan \Jq/p and 6f ^ 0, corresponding to the case of interest, but for now 
let us be general. The solution to (15.81) is then given by 

Ssit) = i6s).cos ^ j + -^=^^(^.).sm 

(5.10) 

where {6s)i and {6s)i are respectively the amplitude and time-derivative of 6s at the onset of 
the transition. From (15. 4p we see that {Ss)i ~ H{6s)i, and therefore {Ss)iT <^ (5s)j. Hence 
the final amplitude of the entropy perturbation is given by 

i6s) f ^ (5s), cos (^VsiOf - Oi)') . (5.11) 



For the case of interest, \6f — 6i\ ^ arctan ^/qjp, which proves our claim that 6s changes by 
a factor of order unity during the transition. 



6 General Remarks on Bouncing Cosmologies 

Next we turn to the all-important issue of reversing from contraction to expansion in a 
non-singular fashion. Now that ( has been shown to be scale-invariant in the contracting 
phase, whatever physics generates a bounce, the usual conservation of C on super-horizon 
scales guarantees that the universe will emerge in the expanding phase endowed with a 
scale-invariant spectrum. 

Unfortunately there is certainly not a wealth of consistent bouncing scenarios to choose 
from. The recently proposed mechanism of Creminelli et al. [H] based on ghost conden- 
sation [T6] is at present the only such mechanism all within the realm of a consistent 
effective field theory and in which all instabilities can be kept under control. There are some 
concerns as to whether something like ghost condensation can be realized in a consistent 
theory of quantum gravity, such as string theory [42j. But for the purpose of this scenario 
we will not be worried about issues of UV completion. The essential point is that it provides 
a consistent and ghost-free effective theory. 

Let us begin by making a few general remarks about obtaining a successful bounce. 
Since H < during a contracting phase, by definition, in order to have a bounce {H = 0) 
we need some epoch during which H > 0. But this is highly non-trivial to achieve in any 
particle physics model. Consider for instance a non-linear sigma model of scalar fields 
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with arbitrary potential [2]: 

= \Gi3 (0") ^0*^0' + ^ (00 • (6-1) 

As long as the metric on moduli space Gij is positive-definite, then a bounce is simply 
impossible, since 

n = -^G^j (0') d(P'W < , (6.2) 

independently of the potential. More generally, for a perfect fiuid with energy density p and 
pressure V, we have 

M|iij = -l(p + P) . (6.3) 

Thus a necessary condition for a bounce is a violation of the NEC: p + V > 0. (In covariant 
form, the NEC is the requirement that T^yU^n'^ > for any null vector [48J.) 

Of course one can get if > by relaxing the positive-definiteness of Gij, i.e. allowing for 
ghost-like modes. But ghosts have disastrous consequences for the viability of the theory. In 
order to regulate the rate of vacuum decay one must invoke explicit Lorentz breaking at some 
low scale ^9]. In any case there is no sense in which a theory with ghosts can be thought as 
an effective theory, since the ghost instability is present all the way to the UV cut-off of the 
theory. More generally, it was argued in that in generic 2-derivative theories, violations of 
the NEC immediately imply the presence of ghosts or tachyons with arbitrarily-large mass. 
The ghost condensate evades this theorem since it relies on higher-derivative kinetic terms, 
while nevertheless defining a consistent effective theory. 

Even if some component violates the NEC without introducing catastrophic instabilities, 
it is still a non-trivial feat to dynamically achieve a bounce. The point is that in order for 
H and then H to reverse sign, this NEC-violating component must come to dominate the 
energy. Since the equation of state is w < —1, by definition, this means that its energy 
density scales like 

Pnec ~ a" , (6.4) 

where a is some positive power. In an expanding universe, this grows while everything 
else redshifts, meaning that this component is bound to dominate the universe. But in 
a contracting universe, of course, precisely the opposite happens. While everything else 
blueshifts, pnec becomes negligible. Put another way, any bouncing solution is unstable to 
the addition of energy into a normal component — in order to achieve a bounce, one must 
ensure that nearly all the energy somehow gets funneled to the NEC-violating fiuid. This is 
a key challenge for any scenario of bouncing cosmology. 
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In the remaining sections of the paper, we show exphcitly how the ghost-condensate 
bounce can be successfully merged with the preceding ekpyrotic phase. The outcome is a 
singularity-free bouncing cosmology in which a scale-invariant spectrum of perturbations is 
generated during the contracting phase and transferred unscathed through the bounce. 



7 Essentials of Ghost Condensation 

Theories of ghost condensation describe a scalar field with higher-derivative kinetic term [TB| 

C = v^M^P(X), (7.1) 

where 

X = -^W)^ (7.2) 

is dimensionless, and M and m are some arbitrary scales to be determined by the fundamen- 
tal theory. (The ghost condensate literature usually defines X = —((90)^/2, with (p having 
dimension of length. Here, however, we stick to the usual mass dimension for scalars, for 
consistency with earlier sections.) Theories of the form (17. ip were first introduced in cos- 
mology as fc-infiation ^50j or /c-essence models [51] to drive accelerated expansion without 
potential energy. 

In a cosmo logical context, the scalar satisfies 
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For generic P{X), this implies the usual redshifting (blueshifting) of scalar kinetic energy as 
the universe expands (contracts). However, if P{X) displays a minimum at some finite X, 
which by rescaling of m can be chosen to lie at Xq = 1/2, then X = 1/2 provides an exact 
solution to (17.31) . This corresponds to the field maintaining constant kinetic energy and thus 
growing linearly in time 

0(t) = -mH . (7.4) 

(Of course this is a solution for any extremum of P{X). However, as we will see shortly, 
fluctuations around a maximum are ghost-like.) The purely derivative nature of the La- 
grangian (17. ip is technically natural if there is some global shift symmetry — + const. 
Moreover, since the above solution is linear in time, corrections involving more than one time 
derivatives on 0, such as (^10)^, vanish identically. Thus, near a minimum of P{X), (17.11) 
provides a consistent effective fleld theory. 
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The full stress tensor of the ghost condensate is derived as usual from (17.11) 

T^u = g^,M^P{X) + Px9^09.0 , (7.5) 
corresponding in the cosmological context to an energy density and pressure given by 

p = M^{2P^xX - P) 

V = M^P{X). (7.6) 

In particular, at an extremum of P the ghost condensate behaves as a fluid with w = — 1, 
mimicking the effect of a cosmological term. 

Ignoring gravity for a moment, small fluctuations in the scalar field, = —m?t + vr, 
around a constant X solution have the quadratic lagrangian 

C cx (Px (Xo) + 2XoPxx (Xo)) TT^ - P^ (Xo) (V7r)2 . (7.7) 

In particular, around an extremum of P(X), perturbations have "right-sign" kinetic term if 
the extremum is a minimum, and are ghostlike for a maximum. Furthermore, the gradient 
term vanishes at that point, meaning that one must keep the leading correction to the P{X) 
lagrangian, say from (□</')^, to obtain 

^grad ~ -^(VV)^ (7.8) 

This implies the dispersion relation uo'^ ~ /M'"^ . (To make contact with the normalization 
convention in [T3], their M is related to our M' by M' = M'^/M.) We postpone the stability 
analysis of small fluctuations, including metric perturbations, to Sec. I7.2[ We first describe 
how NEC violation is achieved in ghost condensation. 



7.1 Violating the NEC 

At the minimum of P{X), the ghost condensate behaves as a cosmological term and as such is 
just at the borderline of violating the NEC. To push it in the w < —1 region, following [13] 
we introduce a potential V{(j)) for the scalar. If this potential is sufficiently flat, then it 
gives a small correction to the scalar equation of motion, allowing us to treat it as a small 
perturbation to (17.41) : (p = —m'^t + TT. To leading order in vr, the energy density and pressure 
become 

p ~ — + v 

V ^ -V, (7.9) 
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where K = P xx(l/2) is the (dimensionless) curvature of the kinetic function about the 
ghost condensate point. Meanwhile we have set P(l/2) = through a constant shift in V. 

And here we discover the culprit for NEC violations in ghost condensate, in the form of 
a term linear in tt in the energy density. Indeed, the H equation — see (16.31) — takes the 
form 

and therefore can have either sign, depending on the dynamics of tt. The latter is determined 
by expanding ^a~^dt (^a^P^x<pj = —y,<j> around the = —iv?t solution: 

This looks like the equation of motion for a scalar field, except that here is to zeroth 
order independent of tt — instead it is a time- dependent source term determined by the 
background (f) = —m^t evolution. 

7.2 Taming the Instabilities 

The ghost condensate model suffers from two types of instabilities: one of the Jeans-type 
which is present even in standard ghost condensation; and a second, gradient-type instability 
which is due to the NEC violating background. We briefly review these in turn, starting from 
the latter. First consider expanding the original P{X) action (17. ip around the NEC- violating 
background 

(j) = -m^t + TTo{t) + Tr{t,x), (7.12) 
where 7io{t) satisfies (17. lip . Using that H n from (I7.10p . the result is [1] 



Thus a non-vanishing H generates a non-zero gradient term for the fluctuations, which evi- 
dently has the wrong sign on the NEC-violating background. Correspondingly, the dispersion 
relation, 

has a gradient instability at long wavelengths. The rate of this instability peaks at k'^ ~ 
HMp^M'"^ /m'^, corresponding to 

2 
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Of course such instabilities are absent if their rate is less than the Hubble rate, [cUgradl ^ 
giving the bound 

The second type of instability arises when including mixing of tt with gravity and is 
present even in standard ghost condensation. One finds a Jeans-type instability with rate [T5| 
[16] 

O^Jeans ~ ^, , (7.17) 

which once again is harmless if less than the Hubble rate: 



<H. (7.18) 



Equations (17.161) and fl7.18p together require that M and M' be chosen such that they 
fit in the range 

Evidently in order to have a parametric window for which this can be satisfied, it must be 
that H ^ H^. If ~ H^, for instance, then (17.191) implies a relation between M, m and 
M'; whereas H ^ implies that one of the instabilities is present. As we approach the 
bounce and \H\ 0, clearly the condition H ^ must break down sooner or later. Since 
the instabilities of interest are of the Jeans-type, however, their effects can be mitigated if 
the entire period of NEC- violation lasts roughly one e-fold, so that 

HAt < 1 , (7.20) 

where H is the Hubble parameter at the onset of NEC-violation, say. 



8 Merging Ekpyrosis and Ghost Condensate 

The next step consists of merging the ekpyrotic phase with the NEC- violating ghost conden- 
sate to reverse contraction to expansion in a non-singular fashion. For simplicity we focus 
first on the case of single-field ekpyrosis — we will discuss the generalization to two fields at 
the end of the section. 

As mentioned in Sec. |6l the existence of an NEC-violating fluid does not guarantee a 
bounce — it must also come to dominate the universe in order to reverse the sign of H. 
But this is a non-trivial feat since its energy density redshifts instead of blueshifting as the 
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universe contracts. Thus, for instance, if the ghost condensate is a spectator field different 
than the scalar relevant for the ekpyrotic phase, one must ensure that somehow the energy 
in the former gets efficiently transfered to the latter, with negligible energy going to other 
light degrees of freedom. This seems to us a highly unnatural possibility. 

Instead if the ekpyrotic and ghost condensate scalars are one of the same field — let us 
denote it again by (p — , then the transfer of energy is of course a non-issue. We will describe 
what this merging entails for the global form of P{X) and and derive consistency 

conditions required for a successful bounce. As we will see, these relations can be satisfied 
for a wide class of kinetic functions and scalar potentials. 

Each subsection deals with different aspects of this merger. For pedagogical purposes, 
we provide at the conclusion of each subsection a short summary of the key results. 

8.1 Global Form of P{X) 

Unlike ghost inflation, where the infiaton is in the ghost condensate phase while driving 
inflation, here it is crucial that the scalar field has approximately standard kinetic term in 
the ekpyrotic phase, i.e. in region a) of Fig. [TJ That is, 

M^P (X) ^ m^X = ~id(j)f (ekpyrotic phase) . (8.1) 

This follows immediately from the fact that V is steep in this region and therefore would 
generically violate the flatness condition of ghost condensate. Furthermore the perturbation 
calculation crucially relies on the kinetic term being canonical. 

Of course in the ghost phase, P{X) must have some minimum where the field can con- 
densate. In its vicinity, we can assume the following quadratic form: 

P{X) ^ :| (^X - (ghost phase) , (8.2) 

so that = —m?t at the minimum while P^xx = K, consistent with the conventions in our 
earlier discussion. 

One might expect that the quasi-linear behavior (18.11) occurs for small X, while the 
minimum in (18.21) is generated when higher- derivative corrections become important at suf- 
ficiently large X. This picture would also be consistent with the fact that X increases in the 
ekpyrotic phase, due to cosmological blueshift. However, a moment's thought reveals that 
such a P{X) would necessarily have a maximum somewhere in between, thereby signaling 
the existence of a real ghost. This is sketched in Fig. [2^. If one is willing to tolerate a real 
ghost, then of course there is no need to invoke all of the ghost condensate technology. 
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Figure 2: Two possible choices for the global behavior of the kinetic function P{X). In case 
a), the linear regime (ekpyrotic phase) lies at smaller X than the ghost condensate point. 
This necessarily implies a maximum for P{X) in between, signaling the presence of a real 
ghost. In case b), the linear regime lies at larger X, thereby avoiding real ghosts. 

So let us instead assume that the quasi-linear regime occurs at larger X than the ghost 
condensate point, as sketched in Fig. [2}d. Thus, at the onset of the ekpyrotic phase, the 
scalar kinetic energy starts out to the right of the minimum, say, and moves away from it. 
Thus the challenge here is to bring the field back to the vicinity of the minimum of P{X) 
at the end of the ekpyrotic phase. As we will see this can be achieved naturally by the form 
of the potential. We will not attempt to motivate this form of P{X) from a UV-complete 
theory. Instead our goal is to provide a specific example of P{X) which, combined with the 
general form of the potential shown in Fig. [1], yields a non-singular alternative cosmology to 
inflation. 

To summarize, we have argued that having a canonical form for the kinetic term during 
the ekpyrotic phase, combined with the no-ghost constraint during the subsequent evolution 
towards the ghost condensate point, requires the kinetic function P{X) to have the form 
shown in Fig. [2)d. 



Since X grows during the contracting phase, it must somehow find its way back to the 
vicinity of the minimum of P{X) at the end of the ekpyrotic phase. The idea is to exploit 
the sharp rise in the potential, shown in region (b) of Fig. [H to greatly reduce the kinetic 
energy in the scalar and bring X near the ghost condensate point. In other words, by 
appropriately choosing the difference in potential energy between the value at the minimum 



8.2 Scalar Potential V{(f)) 
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and that on the plateau, the field will lose just the right amount of kinetic energy to bring 
X near the minimum of P{X). The field evolution in both potential and kinetic function 
space is sketched in parallel in Fig. [31 



P(X) 



PCX) 





We can be more precise about the form of the potential in region c) during the ghost 
condensate phase. Since the latter requires a fiat potential, we can approximate V{(j)) as 
linear in 0: 



^(^)«„A.^1_;3__^J, (8.3) 

where A is some scale, while a and jS are dimensionless. The form of the /5-term is chosen 
for convenience, as will become apparent shortly. 

As we now argue, a and f3 must both be positive. To prove the former one only needs to 
look at the Friedmann equation, which, to leading order in vr, is just 



SH^M^^ = — + — + ah." . (8.4) 



But since tt > during the NEC- violating phase, as seen in (17.101) . it follows that a > 0. 
We can therefore set a = 1 without loss of generality. 

To argue that (3 must be positive as well is equally straightforward. What we want is 
for the NEC to be initially satisfied (tt < 0) and then violated for some time (tt > 0). 
Precisely at the onset of NEC violation, however, corresponding tt = 0, (17.111) reduces to 
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TT = pK^im? / KM'^M-p\. Hence to proceed to the NEC-violating phase, we need /3 > 0. Note 
that a potential of the form (18.31) with a, /? > is indeed consistent with region c) of Fig. [T] 
In this subsection we have shown that the scalar potential in the ghost condensate regime 
must be positive, flat, and have negative slope. The last condition follows from our choice 
that < at the ghost condensate point, to be consistent with the field motion during the 
ekpyrotic phase. 



8.3 Dynamics of the Bounce 

A key approximation that must hold throughout the NEC-violating phase and through the 
bounce for perturbation theory to apply is 

77 < (8.5) 

To see what this entails on the potential and kinetic function, let us combine (I7.10p . (18.31) 
and (18.41) to obtain the master equation 

Ml, (SH^ + 2if) = (^1 + ^ A^ , (8.6) 

valid to leading order in tt. In the last step we have further assumed that f3A'^At/Mp\ < 1 
for the period of interest, as we will prove shortly. Now, since tt H, the condition (18.51) is 
most stringent when H is maximal. Evidently this occurs at the bounce itself, when H = 0, 
and thus H ^ A^/2M|j. Translated in terms of tt, it follows that a necessary condition 
for (18.51) to be satisfied is 

A^<M^K. (8.7) 

But then, looking back at (18.41) . this immediately implies that the Hubble rate is also con- 
strained throughout the ghost condensation phase: 

-<^- 

There is one further condition, namely (17.201) . having to do with the suppression of Jeans- 
like instabilities during the bounce. To see what it implies, let us introduce the dimensionless 
parameters 

in terms of which (18.61) takes the simplified form 

2H + 3H^ = l + pi. (8.10) 
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Figure 4: Evolution of H{t) for (3 = 1. We see that first H vanishes at some time, and 
later on H itself vanishes, marking the reversal from contraction to expansion. Note that, in 
solving flS.lOp . initial conditions were chosen such that t = corresponds to the beginning 
of the ghost condensation phase. 

In terms of these dimensionless variables, the condition fl7.20p is just H^t < 1. Thus, 
provided that 

(3^0{1), (8.11) 

then there is only one time scale in flS.lOp . and thus (17.200 is guaranteed to hold. Figure H] 
shows the result of numerically integrating (18.101) for (3 = 1. The initial conditions are such 
that H <^ at the beginning, ensuring that all instabilities are under control. These 
conditions lead to a NEC- violating phase, followed by a bounce. Since j3 ~ C>{1), the 
duration of the NEC-violating phase is indeed of order one in dimensionless time units, 
satisfying the requirement (17.201) . 

Knowing that (3 ~ C>{1), we can determine the change in (j) during the NEC- violating 
phase: 

2 9 

A(f) = m'^At ~ — ~ TT^Pi ' (8-12) 

where in the third step we have estimated the Hubble parameter during the ghost condensate 
phase by its approximate value at the onset of NEC violation: H ~ A^/Mpi — see (18. 6p . 
We will see further on that m has to be small, therefore (18.121) generically implies that 
moves by a small distance compared to the Planck mass during the ghost condensate phase. 
Incidentally, note from (18.31) that the variation of V is then of order AV ~ A^, and thus 
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V changes by a factor of order one in the process. This justifies the approximation made 
in dHSD- 

To summarize, perturbation theory around the ghost condensate is vahd throughout 
provided that the scale M of the kinetic function is much greater than that of the potential 
as well as the Hubble parameter. Furthermore, in order to keep gradient instabilities under 
control, we have shown that the parameter (3 characterizing the slope of the potential must 
be of order unity. These conditions are easily satisfied for a wide class of potentials and 
kinetic functions. It was explicitly shown that this yields a non-singular bounce. 



8.4 Consistency Relations for a Successful Merger 

Let us now focus on the transition as the field passes through region (b). Since this cor- 
responds to the transition between the ekpyrotic phase and the ghost-condensate bouncing 
phase, we will find consistency relations between the kinetic function P{X) and the scalar 
potential V{(j)). 

As we did when discussing perturbations, we model this transition as instantaneous 
compared to a Hubble time. This is justified since, just before reaching Knin, the field obeys 
the scaling solution and, therefore, is moving rapidly on a Hubble time, as shown by fl2.7l) : 

("3) 

If the rise in V is sufficiently sharp, then (p will indeed reach the plateau in a very short time 
compared to the Hubble time. In this case we can treat H as constant during this transition, 
by energy conservation. Let us denote it by i^min- Then, from (12.71) we obtain 

M^.Hi,^ = --^K,in ~ -pKnin • (8.14) 

Of course this is expected since 0^/2 ^ —V during the ekpyrotic phase. 

Once we reach the plateau, the assumption is that the field is near the ghost condensate 
point. We will see that this constrains Vmin to lie within some parametric window. First, 
to derive an upper bound, we start from the Friedmann equation (18.41) at the onset of the 
ghost phase: 

3Hl,^M^,^ — + A'. (8.15) 

Substituting (18.141) and rearranging, we get 

TT ^ m 7 H 7 . 8.16 

31 



Now a consistency condition for being near the ghost condensate point is that vr <C M^, as 
we recall from fl8.5l) . Since we already know that <^ KM'^ from fl8.7l) . it immediately 
follows that 

Knin < . 8.17 

P 

To derive a lower bound, note that by assumption we must have 0^ ^ during 
the ekpyrotic process, since P{X) has to be approximately linear corresponding to canonical 
kinetic energy. However, we also know that 0^ ~ —2V during the scaling phase, and therefore 

\V\ > (8.18) 

throughout ekpyrosis. Since at this phase V is monotonous we have to require 

iKkl > m\ (8.19) 

where \4k is the value of the potential at the onset of the ekpyrosis. Let us relate \4iin 
and through the number of e-folds M . Since during ekpyrosis the scale factor a stays 
approximately constant, the number of e-folds is 

= (8.20) 

Hence, from the solution (12. 7p it follows that 

iKkl = e-2^|K.in|. (8.21) 

Equations (18.171) . (18.191) and (18.211) together imply that Vmin must lie within the range 

m'e^ « iKninI « (8.22) 

P 

Noting that M is 60 or so, this requires m to be less than M by many orders of magnitude. 
Since these two scales are physically different, the above allowed range can be substantially 
broad. Indeed, from the above considerations, M denotes the cutoff scale by which all higher- 
derivative terms are suppressed, whereas m sets the expectation value of (j) at the minimum 
of P{X). Incidentally, the upper bound is, in fact, conservative since it follows from the 
relatively strict requirement that the field be near the ghost condensate point by the time it 
reaches the plateau. More general transitions will therefore loosen this bound. 

To summarize, in order for the field to land in the vicinity of the ghost condensate point 
at the end of the ekpyrotic phase, the minimum of the potential, V^in, must lie within the 
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above range. The lower bound comes from the requirement that the kinetic energy of is 
approximately canonical and, thus, is large comparing to m. On the other hand it cannot 
be too large. Otherwise by the time (f) reaches the plateau its kinetic function P{X) will not 
be near the condensation point. Since m, M, and K are free parameters at the level of 
model-building, we see that these conditions are satisfied for a wide class of models. 

8.5 Generalization to Two Fields 

We conclude with a few words on the generalization to two fields, which is straightforward. 
Essentially both and ip are assumed to have qualitatively the same P{X) and the same 
qualitative shape for their potentials. Both fields will reach their respective Vmin more or 
less at the same time and will proceed to their respective ghost condensate point. (Of course 
these transitions cannot happen simultaneously since our mechanism for converting entropy 
to adiabatic perturbations relies on one field reaching Vmin before the other. Nevertheless 
this time delay could be negligible compared to a Hubble time, for instance.) Once both 
fields have reaches their ghost condensate points, both act as NEC- violating fluids and drive 
the universe towards a non-singular big bounce. 

9 Discussion of Reheating 

After a non-singular bounce has been successfully completed, the energy in the ghost conden- 
sate must somehow get converted into matter and radiation degrees of freedom in order to 
reheat the universe. We briefly comment on two possible reheating mechanisms, one entirely 
at the level of a 4d effective theory, the other inspired by the original ekpyrotic scenario and 
concepts of heterotic M-theory. 

Let us start with the 4d effective mechanism, which has been mentioned in the ghost 
inflation context [52j. After the bounce, most of the energy of the ghost condensate is stored 
in its potential energy — the kinetic energy is proportional to 2P xX — P, which is small 
near the ghost condensate point. To trigger reheating, one assumes that at some field value 
the scalar potential displays a precipitous drop towards zero, after which it becomes flat 
again. When the field reaches the drop, its kinetic energy will be greatly perturbed away 
from the ghost condensate point. If is coupled to matter fields, this non- adiabatic process 
will excite matter and radiation degrees of freedom. In other words, the potential energy 
difference reheats the universe. At the end of the reheating phase, the scalar field settles 
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back to the ghost condensate point, but with neghgible residual energy. Reheating by this 
mechanism is easily achieved within the context of the new ekpyrotic scenario. 

Much more speculative, the second mechanism is inspired by the original ekpyrotic sce- 
nario. Here, the ekpyrotic scalar field has the geometrical interpretation of the distance 
between a bulk M5 brane and the observable end-of-the-world boundary brane. Reheating 
occurs when the bulk brane inelastically collides and is absorbed by the boundary brane. 
This fusion proceeds through a "small instanton" transition, during which the nature of the 
light degrees of freedom can change; greatly increasing, for example, the number of massless 
scalar fields at the collision. In the context of new ekpyrosis, higher-derivative corrections 
to the M5-brane modulus kinetic term could generate a ghost condensate point, allowing 
for a non-singular bounce before the collision. The subsequent brane collision then excites 
these scalars which, in turn, transfer their energy to matter and radiation and reheat the 
universe. In this context, it is natural to expect that the ghost condensate becomes massive 
and disappears following the small instanton transition. 

10 Conclusion 

In this paper, we have presented a new and fully consistent scenario for the origin of the 
primordial density perturbations. Instead of being generated through a rapid phase of accel- 
erated expansion shortly after the big bang, here the perturbations are generated in a phase 
of slow contraction, long before the big bang. The key breakthrough in this paper is a non- 
singular bouncing cosmology, achieved by successfully merging the ekpyrotic phase with a 
subsequent NEC-violating ghost condensate phase. We have derived the explicit consistency 
relations required for this merger to be successful. These can be fulfilled by a wide class of 
kinetic functions and scalar potentials. The entire cosmological evolution is, therefore, under 
control and can be tracked throughout at the level of a 4d effective theory. 

This framework allows us to settle the controversial issue of the fate of perturbations 
through the bounce. A non-singular bounce allows perturbations to remain in the linear 
regime throughout. More importantly, since the evolution can be described within a 4d ef- 
fective theory, the curvature perturbation, (, is unambiguously conserved and goes through 
the bounce unscathed. To generate a scale-invariant spectrum for ( in the pre big bang 
phase, we have made use of a recently proposed mechanism of entropy perturbation gen- 
eration [25] . This is accomplished by having two ekpyrotic scalar fields rolling down 
their respective negative exponential potentials. We have found that the resulting ampli- 
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tude for ( is remarkably similar to its counterpart in inflation. This confirms earlier claims 
in the ekpyrotic literature that the required level of tuning on the ekpryotic and inflationary 
scalar potentials from COBE normalization is comparable. Additionally, the spectral index 
is slightly red for a wide class of potentials, consistent with recent evidence from WMAP. 

As with inflation, the remaining challenge is to embed this scenario within a UV complete 
theory of quantum gravity, such as string theory. Potentials of the type required in new 
ekpyrosis can be realized in = 1 supergravity, including low energy effective string theories. 
However, the question of whether one can realize ghost condensation in string theory remains 
open. Ultimately the bounce could be generated using an entirely different mechanism. And 
indeed many of the results described here, such as the generation of a scale-invariant ( and 
its propagation through the bounce, would apply equally well. Our motivation in focusing 
on the ghost condensate was to provide as concrete a realization of our scenario as possible. 

At the level of a cosmological scenario, "New Ekpyrotic Cosmology" provides a consis- 
tent alternative paradigm to inflationary cosmology. The two scenarios make distinctive 
predictions for the gravitational wave spectrum: the inflationary spectrum is nearly scale- 
invariant, whereas that of ekpyrotic cosmology is very blue and, therefore, unobservable on 
large scales ^ . Moreover, the generic prediction of the simplest inflationary models is a sig- 
nificant gravity wave amplitude, just below the current sensitivity levels of microwave back- 
ground experiments [6] . Ekpyrosis, on the other hand, predicts an unobservably small ampli- 
tude. Thus the failure to detect B-mode polarization in upcoming experiments would place 
inflation in an uncomfortable corner [B], while lending support to the ekpyrotic paradigm. 
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